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Abstract 

We study the transmission coefficient of a plane wave through a ID finite 
quasi-periodic system - the Frenkel-Kontorova (FK) model - embedding in an 
infinite uniform harmonic chain. By varying the mass of atoms in the infinite 
uniform chain, we obtain the transmission coefficients for all eigenfrequencies. 
The phonon localization of the incommensurated FK chain is also studied in 
terms of the transmission coefficients and the Thouless exponents. Moreover, 
the heat conduction of Rubin- Greer-like model for FK chain at low tempera- 
ture is calculated. It is found that the stationary heat flux J(N) ~ N a , and 

a depends on the strength of the external potential. 
PACS: 05.60.+W, 44.90,+c, 64.70.Rh 

Typeset using REVTeX 



*to appear in Phys. Rev. B 59, Bl No. 13 (1999). 



I. INTRODUCTION 



In recent years, there has been growing interest in studying incommensurate structures 
and commensurate-incommensurate phase transitions in condensed matter physics and dy- 
namical systems^HEl. On the one hand, incommensurate structures appear in many physical 
systems such as quasi-crystals, two-dimensional electron systems, magnetic superlattices, 
charge- density waves, organic conductors, and various atomic monolayers adsorbed on crys- 
talline substrates. One of the simplest prototypes for incommensurate structures is the 
Frenkel-Kontorova (FK) mo delS which describes a one- dimensional (ID) chain of atoms with 
an elastic nearest neighbor interaction, and subjected to an external periodic potential. If 
the mean distance between consecutive atoms is not in rational ratio to the period of the ex- 
ternal potential, the corresponding state is incommensurate. It has been shown by AubryiH! 
in his pioneering work that there exists two different configurations for an incommensurate 
state by changing the strength of the external potential. These two configurations are called 
the sliding phase and pinned phase, and the transformation between these two phases has 
been called the transition by breaking of analyticity. On the other hand, the ground state 
equation of the ID FK model is nothing but the standard map which is widely studied in 
nonlinear dynamical theorylii^Ei. The two incommensurate configurations correspond to the 
invariance circle and the Cantorus, respectively. 

In the past two decades, most of the works about ID FK model have been focused on the 
study of the ground state. Little attention has been paid to the linear excitation or phonon 
modes. As is well known that the excitation of phonons is very important, it determines 
the fundamental properties, such as the wave transmission, heat conduction and other low 
temperature thermodynamics properties, of the underlying material. In turn, the study of 
the transmission coefficients and energy transport gives us insights of the excitation. In this 
paper we shall study systematically these related properties, namely, the wave transmission, 
the heat conduction and the phonon localization in ID FK model. Through this study we 
shall have a clear picture about how the macroscopic phenomena such as the transmission 
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coefficients and heat conduction are related to the localization properties of the phonon 
excitation. Moreover, our study suggests a possible way for experimental observation of the 
phonon gap which characterize the phase transition by breaking analyticity. 

Recently, Burkov et a/.0 have studied the localization properties of the phonon eigen- 
states. They solved the equation of phonon numerically for a FK chain of finite length. It 
is found that the phonon eigenstates are extended and quasiperiodic functions for V < V c , 
whereas for V > V c the eigenstates at band edges of phonon spectrum are more localized 
than that one in the middle of the spectrum, but no exponential localization states have 
been foundEl. In their numerical study, the localization property is inferred by computing 



the participation ratio. Ketoja and SatijaES have studied the eigenfunctions corresponding 
to the minimum frequency w min of the phonon spectra in Cantorus regime by using an exact 
decimation scheme. The phonon eigenstates of the minimum frequency are found to be 
critical for V > V c . 

The method used by Burkov et al is limited to the small system size. In our calculation 
we shall study the wave transmission by making use of the transfer matrix method due to 
the following reason. As is known that there is a nice correspondence between the phonon 
and electron properties. The transmission coefficient has been found to be a very useful 
quantity for the study of localization of electron eigenstatesS. The transfer matrix method 
is a powerful method to study the systems with much larger size. In this paper, we shall 
study systematically the properties of the phonon eigenstates by computing the transmission 
coefficients of a plane wave through a ID FK chain. By using the transmission coefficient, we 
are not restricted to only one frequency, instead we can study the localization properties for 
all frequencies. To quantify the localization, we shall also calculate the Thouless exponents 
in addition to the conventional participation ratio and transmission coefficients. 

Another interesting and fundamental property that is related to the phonon excitation 
is the thermal conductivity. Recently, there has been a renascence interest in heat conduc- 
tion in a variety of ID systemdiiHll, since this problem is essential for our understanding of 
the microscopic origin of the macroscopic irreversibility^. At low temperature, the linear 
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excitations of the underlying system are most important for heat conduction. Rubin and 
Greerf§§ have established the relation between the stationary energy flux and the trans- 
mission coefficients of the phonons. They found that the stationary energy flux approach to 
a finite positive value as the number of atoms goes to infinity for uniform periodic harmonic 
chains, and the thermal conductivity is proportional to N. However, for random mass bi- 
nary harmonic chains, the stationary energy flux is found to be proportional to as the 
number of atoms N tends to infinity, and the thermal conductivity is thus proportional to 
N*. Therefore, how the stationary energy flux and the thermal conductivity depends on 
the particle number N for an incommensurate system is of great interest. This will be also 
investigated in the present paper. 

The paper is organized as the follows. In Sec. |I|, we shall describe the model and 
numerical method for calculating the wave transmission coefficient. The phonon localization 
and heat conduction problems shall be discussed in Sec |T| and Sec. [TV], respectively. A 
brief discussion and conclusion is given in Sec. |V|. 

II. THE MODEL AND NUMERICAL METHOD 



The ID FK model can be described by 



p 2 1 



(1) 



2m 2 

where p n and x n are the momentum and position of the nth atom, respectively. V is the 
strength of the external potential, a is the distance between consecutive atoms without 
external potential. Aubry and LeDaeroni showed that the minimum energy configurations 
are periodic when a/2n is rational (commensurate model) and quasiperiodic when a/2ir is 
irrational (incommensurate model). For an incommensurate FK model, there are qualita- 
tively different ground state configurations separated by the transition by breaking analyticity 
predicted by Aubry. For each irrational a there exists a critical value V c of the external po- 
tential. The V c = 0.9716354 • • llllll corresponds to the most irrational number, golden mean 
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value a/2-K = ^rp^-. Without loss of generality, we restrict ourselves to this particular value 
of a in the numerical calculations throughout the paper. 

The phonon equation describing the physical stability of the atoms in the FK model at 
ground state is 



i) n+ i + i) n -i - [2 + V cos(a£)]^ n = -wVn, 



(2) 



wherex° is the equilibrium position of the nth atom in ground state, ip n is displacement of 
the nth atom from its equilibrium position and uj is the eigenfrequency. 
This equation can be written in the form of transfer matrix 



4>n 



T 



) 



(3) 



with 



T = 

-L r). — 



^ -uj 2 + 2 + V cos(x° n ) -1^ 



V 



o 



(4) 



Thus 



Vn+i 



\ 



N 



N 



with P N = T N T N - 1 ---T 1 . 



(5) 



In order to study the transmission of a plane wave through the ID FK system, we first 
consider a uniform harmonic atom chain with atom mass mo, and the external potential V 
is equal to zero. The eigenstates and eigenfrequencies of this chain are simply 



^ n = A + e iqna + A_e~ iqna and J 1 = ^ sin 2 (iga), 



(6) 



respectively. 

Then we replace the segment between n = 1 and n = N by a finite incommensurate FK 
chain. The atom in the FK chain has mass m. Now we consider an incoming plane wave 
from n = — oo with frequency u = \J~-^ sm(^qa). Thus in the range of n > N + 1 there is 
only outgoing wave. That is 



ih n = Ae iqna + A R e- iqna for n < 

(7) 

%l) n = Be iqna for n>N +1 



From Eqs. (^) and (0), after long calculation, we obtain the transmission coefficient 
B 2 4sin 2 (ga) 
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(8) 



' |-(Piv)iie-^ + (P N ) 21 - (P N ) 12 + (P N )22e iqa \ 2 
where (P^ij, i,j — 1,2 are the elements of matrix P N in Eq. (§). 

If we let the atom mass of the uniform harmonic chain be equal to the atom mass 
in the FK chain, i.e. mo = m, as did in the usual study of electron systems^, then we 
could not obtain the transmission coefficients of the phonons with frequency larger than 
-j=. However, for 7^0, there do exists phonons whose frequencies are larger than 
Therefore, in order to study the transmission coefficients of all eigenfrequencies, we should 
let the mass m of the uniform harmonic chain differ from the mass of the FK chain. In Figs. 
HI (a)-(c), we show the transmission coefficients of the finite FK chains as the functions of 
frequency for different V. The parameter ^- = ls approximated by a converging series 

of truncated fraction: F n /F n+ i (n — 1, 2 ), where {F n } is the Fibonacci sequence. The 

results of Figs. |l](a)-(c) are obtained for iV = Fi 6 = 1597, m = 0.8, and m — 1. In our 
calculations, we first obtained the N atomic positions of the equilibrium ground state of 
the FK chain with free boundary condition, i.e. xq = and xn = lirNa, by the gradient 
method! - !. Evidently, for the plane wave with frequency in the gaps of the phonon spectra, 
the transmission coefficients are zero. From Fig. (|l|), we see that for small V there exists 
a wider frequency range with nonzero transmission coefficients than that for large V. This 
can be understood by the following facts. For V = 0, there exists only one frequency band 
from to = to When V ^ 0, the ground state positions of atoms deviate from that in 
the FK chain without external potential. For small V (< V c ), the ground state positions of 
atoms are periodic or quasiperiodic. Therefore, one band splits into several subbands and 
the bandgaps show up. As V increases, more and more subbands and gaps show up. Also, 
we found that the range of the eigenfrequency becomes smaller as V is increased beyond V c . 
This implies that the eigenfrequencies are attracted somehow for V > V c (see Fig. 1(c)), 



which is similar to the band gaps of the Harper model at the critical pointed. Because the 
eigenstates of the Harper model at critical point are critical, this result is one of the signature 
that the eigenstates of an incommensurate FK chain in Cantorus regime are critical. 

It is well known! that the phase transition by breaking of analyticity is manifested by 
the phonon gap for V > V c . In fact, from Fig. 1, we have seen that as V increased to 1.6 
there is a wide range of frequencies around u = 0, in which the transmission coefficient is 
zero. This is the direct consequence of the appearance of the phonon gap. In order to see 
this transition of an incommensurate FK model, we calculate the transmission coefficient of 
a low frequency wave at different V, which are shown in Fig. (fj). In this figure we plot the 
transmission coefficients of the plane wave with frequency u = 10 -6 as a function of V for the 
FK segments having different lengths. It is obviously seen that there is a sharp decrease after 
V > V c and the decrease becomes sharper for larger N. The corresponding particle number 
are 1597, 2584, and 4181 for curve 1, 2, and 3, respectively. This clearly demonstrates 
the existence of the phonon gap for V > V c . In turn, our results illustrated here suggest 
that the transmission coefficient might be a very good quantity for probing the existence 
of phonon gap in the underlying system. Therefore, measuring the transmission coefficient 
would enable us to detect the phase transition by breaking analyticity experimentally. 



III. PHONON LOCALIZATION 

The above discussed wave transmission is a macroscopic phenomenon. To some extent, 
it reflects the microscopic origin, namely, the phonon excitation in the underlying systems. 
In this section, we would like to study the phonon localization from different approaches. 

Participation ratio -Burkov et. al.@ studied this quantity for the finite FK chain by 
numerically solving Eq. (Q) and computing the participation ratio (PR) 

PR = — lZ ^" ™/ , (9) 

here N is number of atoms. If the eigenfunction is extended, PR tends to be a finite limit 
as N — > oo. If an eigenfunction is exponentially localized, iVxPR will be a finite number 



as the length N goes to infinity. Burkov et al found that for V < V c , all eigenfunctions are 
extended for N up to 987 and for V > V c the PR of states near the band edges decrease with 
increasing N, but the true exponential localization was not observed. Their method has been 
limited by computer memory (RAM). Therefore, in addition to study the participation ratio, 
we shall also study other quantities such as the transmission coefficients and the Thouless 
exponent. Our results given in this section can be regarded as an extension and supplement 
to the previous study by othersS'B. 

Transmission coefficient - The transmission coefficient is also a very good quantity 
reflecting the localization property of eigenstates of electron, and has been used widely to 
study the eigenstates of electron moving in random and aperiodic fields§. If the transmission 
coefficient associated with the eigenstate tends to be a finite limit as iV — > oo, the eigenstate 
is extended. If the eigenstate is exponentially localized, the transmission coefficient will 
decrease exponentially as N increasing. This method has an advantage that it can compute 
transmission coefficients for much larger systems. 

Thouless exponent - Aother important quantity that describes the localization of 
eigenstates is the Thouless exponent. The Thouless exponent for an eigenfunction corre- 
sponding to the eigenfrequency Ui of Eq. (fj) is: 

1 N 

7(^) = iv^£ ln K-^ 2 |- ( 10 ) 

If 7(0^) goes to zero as iV — > oo, then the eigenfunction corresponding to Ui is localized 
exponentially. 

In Figs. (0)-(@), we show the transmission coefficients, the Thouless exponents and 
PR's for the eigenstates of a finite FK chain with iV = 987 for V = 0.4, 1.0 and 1.6, 
respectively. From the numerical results, we find that the transmission coefficients of the 
phonon eigenstates at the band edges (actually, these are quasi band edges because the 
phonon spectra are the Cantor-like sets, the quasi bands consist of many subbands if we 
consider a larger FK chain) are smaller than that in the center of the bands. Because 
the transmission coefficient depends on mo, its absolute value is not meaningful for the 
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study of localization of phonon eigenstates. In order to see if these eigenstates at band 
edges are exponentially localized. We calculate the transmission coefficients for several FK 
chains having different lengths. Fig. |^ shows numerical results for some eigenstates at band 
edges. We can see that these eigenstates are not exponentially localized. This is nicely 
demonstrated by the Thouless exponents shown in Figs. 0(b), 0(b) and |5|(b). There the 
Thouless exponent is about 0.005 ~ 0.02, thus the decay length £ = I/7 ~ 10 2 , which is 
about the order of the size of the FK chains. We also calculated these quantities for larger 
systems (up to 4181, in fact the transfer matrix method allows us to go to even larger size, 
e.g. more than 10,946), but no significant difference has been found. 



Recently, by using the renormalization group transformation, Ketoja and SatijaES have 
studied the eigenf unctions of the minimum frequency for V > V c . They found that phonon 
eigenstates defy localization and remain critical. Furthermore, there exists an infinite se- 
quences of parameter values in the regime of V > V c where the renormalization limit cycle 
degenerates into a trivial fixed point0. 

Here, we would like to get a further picture about the localization for other different sce- 
narios. We shall study not only the transmission coefficients of the degenerated point, but 
also those cases corresponding to the nondegenerated and pseudo-degenerated situation at 
the regime of V > V c . The results are shown in Fig. (|7|). There the transmission coefficients 
as functions of N are drawn for V = 1.75656208382674 • ■ ■, 2.33, and 3.89474285492986 • • ■, 
respectively. These three values correspond to degenerated, nondegenerated and pseudo- 
degenerated cases, respectively. We find that there are not qualitative differences of trans- 
mission coefficients for the degenerated parameter values and other values of V > V c . The 
only minor difference is that the curves for degenerated and pseudo-degenerated points are 
more regular than that for non-degenerated. Another thing worthy to be mentioned is that 
for N fixed, the transmission coefficients decreases as n~ 2 for n < ~ and increases as n 2 for 
n > y [See the insets of Fig. [7| (a)-(c)]. This can be understood by realizing that for small 
part of finite FK chain the positions of atoms looks like irrelative, but not real random, 
the plane wave seems propagate through a pseudo-random medium. But for large part of 
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finite FK chain, there exists a certain correlation between the positions of atoms, and the 
correlation increases as n — > N. The plane wave transports through a correlated media, and 
the transmission coefficients increase as N increased. 

IV. HEAT CONDUCTION 

The properties of the phonon excitations will be also manifested in another macroscopic 
quantity- the heat conduction. It is time now to discuss the transport of energy flux or the 
heat conduction in the finite FK chains. This is a quite interesting problem that attracted 
much attention in recent year jll~§l. In this section, we consider a different model of heat 
conduction, which was discussed by Rubin and Greer§ originally. In this model, the chain of 
N particles (which constitutes the systems) is connected at both ends to semi-infinite chains 
of identical particles. The left- and right-ends are put in thermal equilibrium at temperatures 
Tl and Tr, respectively. The original Rubin-Greer model was for periodic and random mass 
chains with harmonic nearest neighbor interaction. They found that the stationary heat 
flux as a function of N can be expressed in terms of the transmission coefficient: 

J{N) = ^^ j\ 2 N {u)du (11) 

It is shown that J(N) ~ iV for the uniform and periodic chains and J(N) ~ for random 
mass chains. For the FK chains both Tl and Tr must be very low so that the formula (JTTJ) 
can be applied. Fig. @ shows some typical results of j = J(N) / (T L — T R ) for different 
V. In the calculation of transmission coefficients, we let the mass of atoms in the left- and 
right-hand semi-infinite uniform harmonic chains be equal to the mass of atoms in the FK 
chain. We found that the J(N) ~ N for V < V c and J(N) ~ N a for V > V c . For V > V c , 
the a depends on the V. But its range is approximately from 0.83 to 0.87, which is larger 
than 0.5 for random systems. This also implies that the phonon eigenstates of the FK chain 
are extended for V < V c and critical for V > V c . 

It should be pointed out that to express the heat flux in terms of the transmission 
coefficients is valid only at very low temperature, namely the particles oscillate nearby their 
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equilibrium positions. In fact, this is a linearization result of the FK model. In real case, 
the FK model is a nonlinear, thus for the general simulation of heat conduction, one should 
take the approach of molecular dynamics simulationilil. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the transmission coefficient of a plane wave through the 
FK chain by making use of the transfer matrix method. We have been able to calculate 
the transmission coefficients of all phonon frequencies. The localization properties of the 
phonons based on the transmission coefficient and Thouless exponents agree with that by the 
participation ratio. We also studied the Rubin-Greer-like model for the FK chain and find 
that the FK chain likes a periodic chain for V < V c , whereas it looks like a chain somewhat 
between random and periodic for V > V c . Our numerical results confirm that all eigenstates 
are critical in Cantorus regime. This result can be understood as the following. For V < 
V c = 0.9716354 • ■ -, the ground state configuration of atoms is quasiper iodic. The V cos(x^) 
in Eq. (Q) is also quasiperiodic and continuous. It corresponds to the phonon problem of 
Harper equation in the extended regime (V < 2). Therefore, all phonon eigenstates are 
extended. For V > V c , the ground state configuration of atoms is a Cantor-like set. The 
V cos(Xn) takes only some finite values. Consequently, the exponentially localized state does 
not exist, and all eigenstates are criticaS 
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FIGURES 

FIG. 1. The transmission coefficients as functions of frequency u for different V. (a) V = 0.4, 
(b) V = 1, and (c) V = 1.6. The parameters in the calculations are N = 1597, mo = 0.8 and 
|p = ^jp±. The critical value V c = 0.9716543 

FIG. 2. The transmission coefficient of low frequency {uj = 10~ 6 ) wave through the FK chain 
as a function of V. For all three curves mo = 1 and ^ = ^ l ■ The curve 1, 2, and 3 are for 
iV=1597, 2548 and 4181, respectively. 

FIG. 3. The transmission coefficients (a), the Thouless exponents 7 (b) and the participation 
ratio (c) for eigenfunctions of the FK chain at V = 0.4. The particle number of the FK chain is 
N = 987, the mass of the atom is tuq = 0.8, and the winding number is = . 

FIG. 4. Same as Fig. (|) but V = 1.0 

FIG. 5. Same as Fig. (|) but V = 1.6 

FIG. 6. The transmission coefficients of plane waves with frequencies near the band edges 
of phonon spectra as functions of N. The solid circles, squares and diamonds correspond to 
to = 0.369915599017013597, 0.782580479119728611 and 0.950165600964576584 for V = 1.0, re- 
spectively. The solid stars, triangles and crosses correspond to u = 0.907882097032394975, 
1.04046216308972528 and 2.03331741557249979 for V = 1.6, respectively. In all cases m = 0.8. 
The lines are drawn for guiding the eye. 

FIG. 7. The transmission coefficient at the minimum frequency tu m in versus TV* for different V. 
(a) V = 1.75656208382674- • •, (b) V = 2.33, and (c) V = 3.89474285492986 ••• . Inserts are ln(t) 
versus ln(N). 
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FIG. 8. The ratio (j = J(N) / (Tr — Tl)) of the stationary energy flux to temperature difference 
for the FK chain as a function of N for different V. The solid circles, squares, diamonds, and 
triangles are for different values of V =0.4, 1.6, 3.0 and 5.0, respectively. The lines are drawn for 
guiding the eye. 



15 




0.0 0.5 1.0 1.5 2.0 2.5 




CO 



0.0 



-5.0 

ln(t) 

-10.0 



■<■*' 



If 



.< 



(a) 



-15.0 



0.020 

Y 

0.010 



0.000 
0.60 
P %40 
0.20 
0.00 



0.0 



St* 



2 



0.5 



1.0 



A. 



:•• it 

' ''• £ 



-: 
i 

i 

r 



1.5 



* . . 

« : 'J 

"'• i»f 




(b) 



(c) 



•TCj 
CTQ 

H 

o 

CTQ 

CD 
r- K 



2.0 



2.5 



CO 



